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ABSTRACT 

We present high-resolution maps of stars, dust, and molecular gas in a strongly lensed submillimeter galaxy 
(SMG) at z = 3.259. HATLAS Jl 14637.9-001 132 is selected from the if^r^c/z^/-Astrophysical Terahertz Large 
Area Survey (H- ATLAS) as a strong lens candidate mainly based on its unusually high 500 /im flux density 
(~300 mJy). It is the only high-redshift Planck detection in the 130 deg^ H- ATLAS Phase 1 area. Keck 
Adaptive Optics images reveal a quadruply imaged galaxy in the ^-band while the Submillimeter Array and 
the Extended Very Large Array show doubly imaged 880 fim and CO(l^O) sources, indicating differentiated 
distributions of the various components in the galaxy. The stars reside in three major kpc-scale clumps extended 
over '^1.6 kpc, the dust in a compact (~ 1 kpc) region ~3 kpc north of the stars, and the cold molecular gas in an 
extended (~7 kpc) disk ~5 kpc northeast of the stars. The emission from the stars, dust, and gas are magnified 
by ~17, ~8, and ~7 times, respectively, by four lensing galaxies at z ~ 1. Intrinsically, the galaxy is a warm 
(Tdust --40-65 K), hyper-luminous (Ljr - 1.6 x 10^^ L©; SFR - 2000 M© yr"^), gas-rich (Mgas/Mbaryon 
70%), young (Msteiiar/SFR ~ 20 Myr), and short-lived (Mgas/SFR ~ 40 Myr) starburst, without a significant 
active galactic nucleus. With physical properties similar to unlensed z>2 SMGs, HATLAS Jl 14637.9-001 132 
offers a detailed view of a typical SMG through a powerful cosmic microscope. 

Subject headings: galaxies: formation — galaxies: individual (HATLAS Jl 14637.9-001 132) — galaxies: 
interactions 
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1. INTRODUCTION 

Brig ht submillim eter- selected galaxies (SMGs; S^^o > 
5 mJy; iBlain et all ^002) provide a powerful probe into the 
distant Universe. Thanks to the negative ^-correction in the 
Rayleigh- Jeans tail of the dust thermal emission, flux lim- 
ited submillimeter surveys reach an almost uniform integrated 
infrared (IR) luminosity limit across a wide redshift range 
(1 < z < 8) and yield a galaxy population mostly at redshifts 
between 1.7 < z < 2.8 (Chapman et al. 2005). With star for- 
mation rates (SFRs) of 10^ Mq yr"\ the SMGs are the 
most intense star-forming galaxies, despite their inevitably 
short-lived nature (lifetime ^0.1 Gyr). Although such in- 
tense starburst systems are extremely rare in the local Uni- 
verse, SMGs and the Lyman break galaxies may contribute 
equally to the comoving SFR density at z ~ 4 dPaddi et al.l 
[2009 ). In addition to their unique energetics, they also rep- 
resent an important stage in massive galaxy formation. Mul- 
tiple lines of evidence suggest that SMGs are likely the pro- 
genitors of massive elliptical galaxies (e.g., Lillv et al. 1999; 



Swinbank et al. 2006; Arav ena et al.l 12010: ILapi et al.l 1201 II: 
Hickoxet al. 2012), which apparently have formed b ulk of 



their stars rapidly at an early epoch (e.g.. lRenzinill2006|) . 

Our understanding of this important high-redshift galaxy 
population are limited by the sensitivity and spatial resolution 
of current facilities. Gravitational lensing offers an elegant so- 
lution by effectively lifting both limiting factors. Also thanks 
to the negative ^-correction, it is relatively straightforward 
to identify strongly lensed SMGs in large ar ea submillime- 
ter surveys (lBlainlll996l: iNegrello et al.l[20Toh . Extragalactic 
sources with 500 /im flux density 5*500 > 100 mJy are mostly 
strongly lensed or blended SMGs, nearby late-type galax- 
ies, o r radio active galactic nuclei (AGNs) (Negrello et al. 
|2007|) . Fortunately, objects in the last two categories can be 
easily removed using data at other wavelengths, leading to 
an extremely high success rate in identifying strongly lensed 
SMGs with this technique. This simple flux selection has 
produced a few well- st udied strongly le nsed SMGs ("Cos- 
mic Ev elash" z = 2.3: [S w inbank et al. I 2010; LockmanOl 
z = 3. 0: iConlev et al.ll2Qlir [Riechers et al. 201 la; Scott et alj 

[20TTI: ID141 z = 4.2: ^ Coxet al 2011 : Bussmann et al. 2012^; 

HLS J091828.6+514223 z = 5.2: |Combes et al. 2012; and 
HATLAS12-00 z = 3.3, the subject of this paper), which 
were mostly discovered by the Hersche^ Space Observatory 
jPilbratt et al .1120101) . More complex selection processes have 
been proposed (e.g., f Gonzalez-Nuevo et al.l t2012) . which 
would allow selecting hundreds of fainter lensed galaxies with 
Herschel. 

The brightest of t he lensed SMGs could also be detected by 
the Planck mission (Planck Collaboration 2011a). With r^A' 
resolutions, such sources are probably blended with fainter 
sources even in the highest frequency /resolution channels of 
Planck (545 and 857 GHz, or 550 and 350 /im). The Phase- 
I 130 deg^ of the ifgr^' c/z^Z-Astrophysica l Terahertz Large 
Area Survey rH-ATLAS: rEales et al.ll20T0 ) covers 28 Planck 
sources in the Planck Early Releas e Compact Source Cata- 
log (ERCSC; Planck C onaborationl l201 lb). Sixteen of them 
are high Galactic latitude cirrus, ten are low-redshift galax- 
ies, and one is resolved into two similarly bright near by spi- 
rals (NGC 3719 and NGC 3720) (iHerranz et al.ll20T2l) . Only 
one Planck source, PLCKERC857 G270.59+58.52 (^550 = 

Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with important 
participation from NASA. 



1.07 ± 0.12 Jy, ^350 = 1.1 ± 0.2 j£]), is dominated by high- 
redshift candidates. With 18'', 25'', and 36" angular reso- 
lutions at 250, 350 and 500 /im, respectively, Herschel re- 
solves the Planck source into 16 objects within 4.23' ra- 
dius of the Planck position (FWHM = 4.23' at 857 GHz). 
There are 15 faint (^'350 ~ 40 mJy) sources surrounding an 
unusually bright source (HATLAS Jl 14637.9-001 132, here- 
after HATLAS12-00; ^350 = 378 ± 58 mJy). Interestingly, 
the total flux densities of the ^^r^'c/z^Z-detected sources in- 
side the Planck beam account for only 33% and 42% of 
the Planck flux densities at 545 and 857 GHz, respec- 
tively, suggesting that the Planck measurements are overes- 
timated bec ause of blending with an over-density of sources 
(INegrello et al. 2005) or positive noise spikes (i.e., "Edding- 
ton Bias";lddington 1913). HATLAS 12-00 peaks at 350 /im 
in flux density ("350 /im peaker"), implying a high photomet- 
ric redshift given typical dust temperatures. Subsequent de- 
tections of multiple carbon-monoxide (CO) lines from this un- 
usually bright object determined a redshift of zco = 3.2592 
0.0010 (Zspectrometer, CARMA, Z- Spec: . Harris et al.ll2012L 
Van der Werf et al. in prep; Riechers et al. in prep). 

The high 500 /im flux density (^500 = 298 ± 24 mJy), 
in combination with the confirmed high redshift, makes 
HATLAS 12-00 an excellent strong lens candidate. In this 
paper we present a detailed multi-wavelength analysis of 
this P/a^c^-associated SMG. We describe our high-resolution 
Keck adaptive optics imaging, Submillimeter Array (SMA) 
and Expanded Very Large Array (EVLA) interferometric ob- 
servations, and the panchromatic photometry in § [2l We then 
perform a joint strong lens modeling at rest- frame 0.5 /im, 
200 /im, and CO(l^O) in § [51 In § H we derive the in- 
trinsic physical properties of the SMG from its spectral en- 
ergy distribution (SED). We conclude discussing the physical 
properties of HATLAS 12-00 in the context of unlensed z > 2 
SMGs (§[5]). Throughout we adopt a ACDM cosmology with 
l^m = 0.3, VLa = 0.7 and Hq = 70 km s'^ Mpc-\ 

2. OBSERVATIONS 
2.1. Keck Adaptive Optics Imaging 

We obtained a 3,440-s ^s-band (hereafter K) image on 201 1 
April 13 (UT) and a 2,100-s /-band image on 2011 June 30 
(UT) with the K eck II laser guide- star adaptive-optics system 
(LGSAO; .Wizinowich et al.l 120061) . An R = 15.8 magnitude 
star 48" SW of HATLAS 12-00 served as the tip-tilt refer- 
ence star. The estimated Strehl ratios at the source position 
are ~23% and 5% in ^ and /-band, respectively. We used 
the NIRC2 camera at 0."04 pixel"^ scale for both filters (40" 
field), and dithered with 2-3" steps. The outside seeing at 
0.5 /im was ^0.4" and 0.5" during the K and /-band imag- 
ing, respectivel}|3. 

We used our IDL (Interactive Data Language) programs to 
reduce the images. After dark subtraction and flat-fielding, 
sky background and object masks are updated iteratively. For 
each frame, after subtracting a scaled median sky, the residual 
background is removed with B- spline models. In the last iter- 
ation, we discard the three frames of the poorest image quality 
and correct the NIRC2 geometric distortion using the solu- 
tion of R B. CamerorE3 before combining the aligned frames. 
The resolution of the final K and /-band images are 0."16 

The source is only detected at 545 and 857 GHz by Planck. 
http://kiloaoloa.soest.hawaii.edu/current/seeing/ 

http://www2.keck.hawaii.edu/inst/nirc2/forReDoc/post_observing/dewarp/ 



A Strongly Lensed SMG at z 3 3 



I 40 kpc (z=1) I 



Peak Positions 



SMA'880um 



V (10^ km/s) \ 
-6 -4 -2 2 4^ 



-^EVL<§P(l"-Ot- 



1 




Fig. 1. — High-resolution images of HATLAS 12-00. All images are aligned and the tickmarks are spaced at intervals of I" . Green crosses mark the two 
components seen in the EVLA image, a: Keck K-bmd image painted with a pseudo-colormap from Keck K (Red), / (Green), and ACAM optical (Blue) images. 
Lensing galaxies and the PSF star are labelled. The scale bar indicates 5^^ or 40 kpc at the lens redshift. The inset shows the lens -subtracted K-bmd image 
overlaid with the peak positions for lens modeling (§ 13. IV For clarity, the positional errors, as indicated by the ellipses, are enlarged by a factor of four. The 
colors distinguish images from the three clumps in the source plane, b: SMA 880 iim image. Contours are drawn at -2,-1,-1-1,-1-2, and +4a, where a is the r.m.s. 
noise (3 mJy beam~^). c: EVLA CO(l— )^0) image. Contours are drawn at -1,-1-2,-1-4, and -i-8cr, where a is the r.m.s. noise (27 /iJy beam~^). The inset shows the 
CO spectrum from the same data cube, along with a Gaussian fit (red). In b and c, the ellipse to the lower right shows the beam. 



and 0.''27 in FWHM, respectively. We measure the FWHMs 
from the most compact source in the field located 10'' SE 
of HATLAS12-00 (labeled "PSF" in Fig. [Tfc); we also use 
this object as the PSF in the lens modeling (§ 13.21) . The im- 
ages are fl ux calibrated a gainst UKIRT Infrared Sky Survey 
(UKIDSS: lLawrence et al.l2Q Q7) and reach depths of ^ = 25.6 
and / = 25.0 AB for a 5cr detection with O'.'l and 0'/2 radius 
apertureffl respectively. 

2.2. William Herschel Telescope Imaging 

Limited by the small field of NIRC2, a deep wide-field im- 
age is required for astrometry calibration. Optical imaging 
was obtained with the high-throughput auxiliary-port camera 
(ACAM) mounted at a folde d-Casseg rain focus of the 4.2-m 
William Herschel Telescope (iBenn et al. 2008) on 201 1 April 
26 (UT). We obtained four images of 200 s on a ~2' field 
centered on HATLAS 12-00, without any filter. The seeing 
was ~0.''9. The images were reduced and combined follow- 
ing standard techniques in IRAlQ. No accurate photometric 
calibration is possible because we did not use any broad-band 
filter. But by comparing sources extracted from the ACAM 
image and the SDSS /-band catalog in the same field, we find 
that our image reaches an equivalent /-band 5 -a depth of 24.6 
AB, or 2.3 magnitudes deeper than the SDSS. 

We solve the astrometry of the ACAM image using 
the on-sky positio ns of Sloan Digital Sky Survey (SDSS; 
lAihara et al.(l2011h DR8 sources inside the field. We use the 
astrometry routines in Marc Buie's IDL librar}0 to correct 
for offsets, rotation, and distortions, with four terms (con- 
stant, X, y, and R = Vx^-\-Y^). Sources that appear blended 
in the SDSS catalog are excluded. With 35 SDSS sources, we 
measure la dispersions of SRA = 0.''13 and (5Dec = 0.''14 be- 
tween the astrometry calibrated ACAM image and the SDSS. 
Finally, we use the same routines to match the NIRC2 images 
to the ACAM image with 13 well-detected sources inside the 
40'' NIRC2 field of view. The corrected NIRC2 images show 

Different aperture sizes were chosen here because of the different reso- 
lutions. 

http://iraf.noao.edu/ 
http://www.boulder.swri.edu/~buie/idl/ 



Icr dispersions of 5RA = 0."04 and (5Dec = 0."05. 

2.3. SMA 880 jj^m Imaging 

We obtained SMA interferometric imaging of 
HATLAS12-00 at 880 /im on 2011 May 2 with 4 hrs 
on- source integration time. The SMA single-polarization 
receivers were tuned with the primary local oscillator at 
339.58 GHz and an intermediate frequency coverage of 4-8 
GHz, providing a total of 8 GHz bandwidth (considering 
both sidebands). The SMA was in its compact configuration, 
providing a FWHM resolution of 2."07 x l."87 at PA = 
-23.6° (Position Angle, E of N) with a primary beam of 
-37". The quasars 1229+020 and 1058+015 were used for 
time- variable gain (amplitude and phase) calibration. The 
blazar 3C 279 was used as the primary bandpass calibrator, 
while Titan was used as the absolute flux calibrator. 

The SMA image has an r.m.s. noise level of 3 mJy beam~^ 
and a pixel size of 0."1. The Herschel source is resolved 
into two components (Fig. [TJ: the southern one (RA = 
11:46:37.997, Dec = -00:11:33.53 [J2000]) is slightly re- 
solved and has an integrated flux density of 43 ± 6 mJy, 
and the northern component (RA = 11:46:37.951, Dec = 
-00:11:30.32) is point-like and have a peak flux density of 
27 ± 3 mJy. We therefore adopt a total flux of 70 ± 10 mJy, 
taking into account the 10% flux calibration uncertainty. 

2.4. EVLA C0( 1^0) Imaging 

We exploited the recent upgrade to the Nat ional Radio 
Astro nomy Observator}|3's Very Large Array (iPerley et al.l 
1201 1|) . which includes the provision of Ka-band re- 
ceivers (26.5-40 GHz), to observe the redshifted CO(l^O) 
emission from HATLAS 12-00 at 27 06532 GHz (i^rest = 
115.27120GHz: lMorton & Noreaul[T994l) . 

Observations were carried out dynamically during excel- 
lent weather conditions on 2012 January 6 and 8. During this 
Open Shared Risk Observing period the available bandwidth 
from the new Wideband Interferometric Digial ARchitecture 

NRAO is operated by Associated Universities Inc., under a cooperative 
agreement with the National Science Foundation. 
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(WIDAR) correlator consisted of two independently tunable 
output pairs of eight sub-bands each, with 64 x 2-MHz full- 
polarisation channels per sub-band, giving a total bandwidth 
of 2,048 MHz. At the redshift of HATLAS 12-00, however, 
the CO(l^O) line could be reached by only the BD output 
pair, giving ~ 11,350 km s~^ coverage and ~ 22 km s~^ reso- 
lution. We offset our tuning by 64MHz to avoid noisier edge 
channels. The 8 sub-bands of output pair AC were tuned to 
32.5 GHz. 

The bright compact calibration source, Jl 150-0023 were 
observed every few minutes to determine accurate com- 
plex gain solutions and bandpass corrections. 3C286 (S = 
2.1666 Jy at 27.06 GHz) was also observed to set the abso- 
lute flux scale, and the pointing accuracy was checked locally 
every hour. In total, around 2hr of data were obtained for 
HATLAS 12-00, with -1 hr of calibration. 

The data were reduced using AIVS (3ip EC12) follow- 
ing the procedures described by llvison et al.l (|201 1). though 
with a number of important changes: data were loaded us- 
ing bdf2aips and fring was used to optimize the delays, 
based on 1 min of data for 3C 286. The base bands were knit- 
ted together using the NOIFS task, yielding uv datasets with 
512 X 2-MHz channels, which we then added together us- 
ing the task DBCON. Finally, the channels were imaged over 
a 512 X 512 X 0.3'^ field, with natural weighting (ROBUST = 
5), to form a 512^ cube centered on HATLAS12-00. Inte- 
grating over those 55 channels found to contain line emission 
(so a FWZI of ~ 1, 200km s~^) yielded an r.m.s. noise level of 
27 /iJy beam~^. 

The CLEANed and velocity-integrated CO map is shown 
in Fig. [It. The beam is 2f'5x2f'2 at PA = 85°. Simi- 
lar to the SMA, the map resolves two components separated 
by ~5'^ The CO lines extracted from the two components 
show the same redshift and line profile, further confirming 
that they are lensed images of a single source. The best-fit 
Gaussian to the area-integrated spectrum gives a line width 
of AVfwhm = 585 ± 55 km s"^ and a line flux of ^coAV = 
1.52±0.20Jy km s~^ . In comparison , the CO(l^O) measure- 
ments reported by lHarris et al.l (12012) using Zpectrometer on 
the Green Bank Telescope are: AVfwhm = 680 ± 80 km s~^ 
and ^coAV = 1.18 ± 0.26 Jy km s-^corrected for the 20% 
difference in the absolute flux density of 3C 286). Here- 
after, we use the weighted mean of the two measurements, 
Sco^y = 1.40 ± 0.22 Jy km s~^ to derive the molecular gas 
mass. 

2.5. Panchromatic Photometry 

Photometry of HATLAS 12-00 were obtained from the 
SDSS (w,g,r,/, and z), the UKIDSS (Y,J,H, and K), the 
Wide-Fiel d Infrared Survey Explorer (WISE, 3.6 and 4.6 /im; 
IWright et al. 2010), the Herschel/FACS (100 and 160 /im; 
Program ID: 0T1_RIVIS0N_1; (Ib^etal. 2010), th e Her- 
schel/SFlRE (250 , 350, and 500 /im; fPascale etlH I20TTI: 
iRigbv etal.l 1201 Ih . the Large APEX BOlometer CAmera 
(LABOCA, 870 /im:'Siring o etani2009h . the SMA (880 /im), 
the Max-P lanck Millimetre Bolometer (MAMBO, 1.2 mm; 
[Kreysa et al. 1999), the Combined Array for Research in 
Millimeter-wave Astronomy (CARMA, 2792 and 3722 /im; 
iBocketal. 2006), and the VLA FIRST survey (21 cm; 
iBecker et al 1995) . 

We obtained imaging at 870 /im with the LABOCA 
bolometer array at the Atacama Pathfinder Experiment 
(APEX) telescope in November 2011 (Clements et al. in 



TABLE 1 
Photometry 



Instrument 


Band 


A 


F^(G1+G2) 


F^(SMG) 






(/im) 


(/iJy) 


(mJy) 


SDSS 


u 

8 


0.36 
0.47 


0.1 zb 0.1 
0.4 ±0.4 






r 
i 


0.62 
0.75 


3.3 ± 1.2 
8.6 zb 2.0 






z 


0.89 


23 zb 8 






Y 


1.03 


35 zb 6 






J 


1.25 


67 zb 8 


0.0017 ± 0.0003 




H 


1.63 


68 zb 10 






K 


2.20 


139 ± 10 


0.0123 ±0.0009 


WISE 


w\ 


3.35 


205 ± 20 


0.037 ±0.020 




w2 


4.60 


242 ±53 


< 0.117 




w3 


11.56 


735 ± 35 


< 0.702 




w4 


22.09 


< 3460 




PACS 


green 


100 




25±6 




red 


160 




138±21 


SPIRE 


blue 


250 




323 ±24 




green 


350 




378 ±28 




red 


500 




298 ±24 


LABOCA 




870 




103 ±19 


SMA 




880 




70 ±10 


MAMBO 




1200 




38±6 


CARMA 




2792 
3722 




1.4±0.5 
<2.0 


VLA 




214000 




1.2±0.4 



prep). LABOCA observed a 11.4' diameter field with a reso- 
lution of FWHM = 18.^^6. The observations have a total inte- 
gration time of '^30 hr reaching a la sensitivity of ^2 mJy. 

We obtained 1.2 mm imaging with MAMBO at the IRAM 
30-m telescope (FWHM ~ 10.^^7) in January and February 
201 1 (Dannerbauer et al. in prep.). Observing time in the on- 
off mode is 24 minutes, achieving a la sensitivity of ~1 mJy. 

We obtained continuum observations at 81.2 and 
108.2 GHz (3722 and 2792 /im; covering rest- frame 
CO[3^2] and CO[4^3] lines) on 2011 Mar 18 and Sep 
1 as part of our CO follow-up campaign of bright, lensed 
H-ATLAS SMGs with CARMA in D array (Riechers et al. 
in prep). Observations were carried out for 0.9 and 1.4 hr 
on source, respectively, using the 3 mm receivers and a 
bandwidth of 3.7 GHz per sideband. HATLAS 12-00 is 
unresolved in these observations, with angular resolutions of 
6'/^ X 5'/0 and 6'/0 x 3f'8 at 81.2 and 108.2 GHz, respectively 
(restored with natural baseline weighting). 

Table [T] lists the photometry. We have included in the er- 
rors the absolute flux calibration uncertainties (3% for WISE, 
3-5% for PACS, 7% for SPIRE, 10% for SMA, and 15% for 
LABOCA, MAMBO, and CARMA). 

3. LENS MODELING 

Because the LGSAO image has the highest spatial resolu- 
tion, we use it to find the best-fit lens model. We initially use 
the peak positions of th e mu ltiply-imaged source to constrain 
the lensing potentials (§ 13.11) , then we exploit the ^-band light 
distribution in the image plane to quantify the morphologies 
of the source as well as refining the lensing potentials (§ 13. 2b . 
Finally, we use the best-fit lensing potentials and the SMA 
and EVLA images to constrain the sizes and locations of the 
dust and molecular gas in the source plane (§ l3.3l & [34b . 

3.1. K-band Peak Positions 

We use LENSTOOL (K neib et al.ll 19961: 1 Julio et al.ll2007h to 
find the best-fit parameters and their errors from the peak po- 



A Strongly Lensed SMG at z 3 



5 



sitions. LENSTOOL implements a Bayesian Markov chain 
Monte-Carlo sampler to derive the posterior distribution of 
each parameter and an estimate of the evidence for the model. 

The tensing system is mainly made of two red filaments 
that are ~3.''5 apart (Fig. [T]^). The outward curved shape 
of the northern arc can be explained if the source is intrin- 
sically curved. Hence we split each of the two arcs into three 
parts and build a simple lens model by putting two deflectors 
centered on Gl and G2. We find that the predicted counter- 
images can explain the additional features close to Gl and 
G2. Guided by the predicted counter-images, we define three 
systems of lensed images (Fig.[T]2 inset). The 11 peak posi- 
tions in 3 separate systems provide a total of 16 constraints 
(11 X 2-3 X 2), allowing us to include shear from nearby 
galaxies G3 and G4. 

For the tensing galaxies, we find photometric redshifts of 
ZG1+G2 = 1.06 ±0.16 and z g4 = 0.80 ±0.28 with the public 
photo-z code EAZY (Bramme r et al.ii200 8). We obtain the 
nine-band photometry from the SDSS (w,g,r,/, and z) and 
the UKIDSS (Y,J,H, and K) surveys. At these wavelengths, 
the flux from the lensed galaxy is negligible (§ H]). G3 is 
undetected in SDSS but shows similar color as Gl and G2. 
Hence in the lens modeling we assume all four galaxies are at 
z=1.06. 

For the lens model, we assume that the dark-matter plus 
baryonic mass profiles of the foreground lens galaxies Gl to 
G4 can be described as singular isothermal ellipsoids (SIEs; 
iKormann et al.l [T994h . The SIE profile is parametrized by 
the velocity dispersion (a), the position (x, y), the axis ra- 
tio {q = b/a), and the PA (6>, E of N). We fix the positions 
to the centers of the galaxies. For G3 and G4, we further 
fix their q and to those from the light distribution, because 
they are not well constrained by the peak positions and there 
are significant correlations between the PA and ellipticity of 
the light and of the mass distribution (e.g. Sluse et al. 2011). 
Therefore, we have a total of eight free parameters. We find a 
best-fit with = 7.9 for dof = 8 (degrees of freedom) and an 
average positional error of 0.''04 (~1 pixel). The parameters 
and their errors are summarized in Table [21 We also list the 
mass enclosed by the critical curve for each SIE, 



TABLE 2 

Derived Properties from Lens Modeling 



Me 



4^^ DlDls (j' 
G Ds c2 



(1) 



where Di, Ds, and Dis are the angular diameter distances to 
the lens, to the source, and between the lens and the source, 
respectively. The radius of the area enclosed by the critical 
curve can be approximated by the circularized Einstein radius: 



b = 4.5( 



200 km s- 



2q_ 
\+q' 



kpc. 



(2) 



In the errors of masses and velocity dispersions, we have in- 
cluded the 1-cr uncertainty of the photometric redshift. 

The nominal model described above is an optimal descrip- 
tion of the tensing system because of the following: 

1. Adding ^'s and 6>'s of G3 and G4 as free parameters 
does not substantially improve the fit: the Bayesian evi- 
dence increases by only A \n{E) = 0.5E3 and the reduced 
actually increases from 1.0 to 1.6 as a result of the 
decreased degree of freedom. 

The improvement of a model is substantial if 1 < Aln(£^) < 2.5, strong 
if 2.5 < A\n(E) < 5, and decisive if Aln(£) > 5 (Jeffrevs 1961). 



Object Quantity 



Value 

Peak Position^ Surface Brightness^ 



Unit 



Gl 


q 




0.46 ±0.10 
-59 ±10 


0.52 ±0.04 
-61 ±2 


deg 


G2 


a 

q 
e 


209 ±24 
1.6+0-9 X lOii 
0.63 ±0.06 
59±7 


212 ±24 
1.7+0-9 X lOii 
0.69 ±0.03 
67±2 


km s~^ 
M© 

deg 




G 


240 ±28 

2.r\i X loii 


240 ±27 
2.7!}-5 X 10^1 


km s~^ 
M© 


G3 


cr 


243 ±51 
2.9+11 X 10^^ 


242 ±28 
2.8!}-^ X 10^1 


km s ^ 

M© 


G4 


G 


168 ±40 
6.6!8-9 X lOio 


165 ±20 
6A+1-1 X 10^0 


km s~^ 
M© 


SMG 


/i(880/im) 
/i(CO) 




16.7 ±0.8 
7.6±1.5 
6.9±1.6 





Not e. — a - Best-fit parameters from the K-bmd peak positions 
(§ 13. U . b - Be st-fit parameters from the K-bmd surface brightness dis- 
tribution (§ 13.2b. c - Total mass enclosed by the critical curve of each SIE 
potential (Eqlflj&EJ. 



2. Excluding the potentials of G3 and/or G4 does degrade 
the fit significantly. The Bayesian evidence decreases 
by Aln(£^) = 2.5 and 30, and the reduced increases 
from 1.0 to 1.8 and 6.5, when we exclude G4 and both 
G3 and G4, respectively. 

3. Including a group-scale potential with a PIEMD pro- 
file (Pseudo-Isothermal Elliptic Mass Distribution; 
iKassiola & Kovnei]ll993 ) does not improve the fit. For 
the PIEMD profile, we adopt a cut-off radius of 500 kpc 
but allow the position, ellipticity, PA, core radius, and 
velocity dispersion to vary. To limit the number of free 
parameters, we fix the ^'s and PA's of the SIEs to those 
measured from the tensing galaxies but allow their ve- 
locity dispersions to vary. So we have a total of 10 free 
parameters. The reduced x^ of the best-fit is 1.9, much 
higher than that of the nominal model. The Bayesian 
evidence also decreases by Aln(£^) = 6.9 when com- 
pared with the nominal model. 



3.2. K-band Source 

Although the peak positions can constrain the deflectors 
through ray tracing, they cannot provide an accurate estimate 
of the magnification factor because the source-plane light dis- 
tribution is not taken into account. Because we want to esti- 
mate the intrinsic properties of the lensed galaxy, we are inter- 
ested in the luminosity weighted magnification factor, which 
depends on the source morphology because the magnification 
factor is different at each source plane position. In this section 
we model the morphology of the lensed galaxy and refining 
the tensing potentials simultaneously with the ^-band image. 
The PSF is derived from the most compact source in the field, 
which is 10'' SE of HATLAS12-00 (Fig.[B/). 

Following § 13.11 we assume that the source consists of three 
clumps, each described as a Sersic profile. Again we use SIE 
profiles for the tensing potentials. So we have a total of 29 pa- 
rameters: seven parameters for each Sersic profile, and eight 
parameters for the SIE potentials. Our fitting procedure is 
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as follows. For an initial s et of parameters describing the 
source and the lenses from § 13.11 we use lenstool to gen- 
erate lensed images of the source, which is then convolved 
with the PSF and compared with the observed image. We 
limit the comparison in a ?>'.'?> x A'.' 5 (83 x 113 pixels) rect- 
angular region that encloses the lensing features. This pro- 
cess is iterated with AMOEBA_SA to find the parameters 
that minimize the residual between the observation and the 
model. AMOEBA_SA is based o n the IDL multidi mensional 
minimization routine AMOEBA (iPress et al.ll 19921) with sim- 
ulated annealing added by E. Rosolowsky. We allow a max- 
imum of 1000 iterations in each call of AMOEBA_SA. For 
the simulated annealing, we adopt an initial "temperature" 
of 100 and decrease it by 20% in each subsequent call to 
AMOEBA_SA. A good fit with a reduced around unity is 
normally found after a few calls to AMOEBA_SA. For each 
iteration, we compute the total luminosity-weighted amplifi- 
cation factor {(Ik) by summing the pixel values in the image 
and the source planes with apertures matched by inverting the 
image plane aperture to the source plane. The l-cr confidence 
interval of (Ik is found with X^(M)~Xmin — ^- ^^^^ 
compute the values on the residual image binned by 4-pixel 
boxes (FWHM = = A pixel), so that the noise becomes 
uncorrected between pixels; or equivalently, one could divide 
the x^ values from the original residual images by a factor of 
16. We find the luminosity-weighted magnification to be (Ik = 
16.7 ± 0.8. The best-fit parameters for the deflectors are listed 
in Tabled The results are very similar to those from fitting the 
peak positions, although the errors are smaller because the en- 
tire image provides more information than the peak positions 
alone. 

Figure [21 shows the best-fit model. The lensed galaxy has a 
curved morphology, causing the northern arc bending in the 
opposite direction of the deflectors. In the source plane, the 
three clumps extend over only 0'/2l or 1.6 kpc, and their ef- 
fective radii are 0f'21 ±0f'04 (1.5 ±0.3 kpc), Of'085 ±0f'013 
(0.6±0.1 kpc), and 0f'll±0f'05 (0.8±0.4kpc) from W to 
E. 

The nature of the feature NNE of G2 in the residual im- 
age is unclear, but it is unlikely to be at the same redshift 
as the lensed galaxy: tracing its position to the source plane 
and imaging it back predicts an unobserved equally bright 
counter-image O'/^S of the southern arc. We believe this fea- 
ture is part of the galaxy G2. 

3.3. 880 fim Source 

Precise astrometry calibration is crucial for a joint anal- 
ysis of images from different wavelengths. Because the 
only K source detected by the SMA is HATLAS 12-00, we 
have to estimate the astrometry offset between the two im- 
ages in a statistical way. Because the Keck image is tied 
to the SDSS astrometry and the SMA image is tied to the 
radio re ference frame, we cross-correlate the VLA FIRST 
catalog (iBecker et al.lll995h and the SDSS catalog within 1° 
of HATLAS 12-00 and compute the optical-radio separation. 
Ninety-four radio sources have optical counterparts within 3". 
We then fit an elliptical Gaussian to the two-dimensional dis- 
tribution in ARA = -(apiRST - <^sdss) and ADec = (5first - 
^^SDSS- The systematic offset from the peak position of the 
Gaussian is consistent with zero (ARA = -0.''07, ADec = 
+0.''08). The best-fit Gaussian has a's of 0.''40 and 0.''30, 
and a PA of 119° for the major axis. Therefore, the la el- 
lipse of the astrometry offset has major/minor semi-axes of 



0.^^61/0 .^^45. Our result is consistent with that of llvezic et al.l 
(|2002|) . who found a ~0.''l systematic offset and a Icr error 
circle of 0.''47 in radius between FIRST and SDSS astrome- 
try. 

We can constrain the astrometry offset further through lens 
mo deling. As demonstrated by Kochanek & Narayan (19921) 
and lWucknitzl (|2004l) . interferometric data are most naturally 
modeled with the i/v-plane visibilities, because it avoids beam 
deconvolution and naturally handles correlated noise. Here, 
however, we opt to model the CLEANed map directly, because 
(1) the images are essentially unresolved in the SMA map, 
and (2) we al ready have a good lens model from the ^-band 
image (§ 13.21) . Because of the limited spatial resolution of the 
SMA 880 /im image, the two centroid positions do not offer 
enough information to constrain the lens model. Hence, for 
the deflectors, we fix all the parameters to the best-fit values 
from § 13. 2[ for the source, we assume a circular Gaussian pro- 
file with variable position and size. We shift the SMA image 
relative to the K model on a 1" ^1" grid with O.'^l steps. At 
each offset position, we find the best-fit model using the same 
fitting procedure as in § 13.21 The modeling is performed on a 
51 X 68-pixel (5.''1 x 6.''8) region enclosing the SMA sources. 
Figure [3]2 shows a map of the minimum x^ values at each off- 
set position. The global best-fit, with reduced of unity, is 
reached when we shift the SMA image 0.''6 E of the K image. 
The middle panels of Fig. [2l show this global best-fit model. 

The noise of the SMA map is Gaussian but is highly cor- 
related. We compute the r.m.s. noise of the SMA map af- 
ter binning it by boxes of ri^ pixels. We find that the noise 
starts to decrease as 1 for ^ > 20 pixels (FWHM 1" = 
20 pixel), indicating that the noise becomes uncorrected on 
20-pixel scales. Therefore, we divide the x^ values from the 
the residual images by a factor of 400, which is equivalent to 
computing x^ from residual images binned by 20-pixel boxes. 

In combination with the la error ellipse from FIRST- 
SDSS cross-correlation, we determine that the astrometry off- 
set between 880 /im and K images is ARA = -0.''5 ± O'/l 
and ADec = O'/O ± 0'/2; i.e., the overlapping region between 
the ellipse and the la contour of the x^ niap. Collecting 
all of the solutions in this permitted offset region satisfying 
X^(/i)-Xmin ^ 1. we estimate a luminosity- weighted 880 /im 
magnification of /iggo = 7.6 ± 1.5, and an 880 /im source size 
of FWHM = Of'15!g 14 = 1.2+10 kpc. Because we have fixed 
the deflectors with the best-fit parameters from ^-band, the 
errors here do not include the uncertainties of the deflectors. 
Higher resolution far-IR images are required to constrain the 
deflectors and the source simultaneously. 

Dust emitting regions are often spatially offset from the 
UV/optical emitting regions in SMGs (Tacconi et al. 2008; 
ICarilli et aLllMol: llvison etaHIMlOa: Riechers et a l. 2010). 
This is clearly the case for HATLAS 12-00, which shows dis- 
tinctly different morphologies at ^-band and 880 /im, even 
after convolving the ^-band image with the SMA beam (com- 
pare Fig. [2]: & /). From the lens model, we estimate a source- 
plane separation between the 880 /im source and the central 
K clump of 0.^^41 ± O'/Ol or 3. 1 ± 0.5 kpc (Fig. [2^). 

If we assume zero astrometry offset between SMA and 
Keck, then we obtain a model that poorly fits the observation 
(Ax^ ~ 4; Fig. The lens model gives a slightly larger 
magnification (/iggo = 8.4 ± 1.6) and doubles the source size 
(FWHM = 2.5+Q3 kpc). However, the source-plane separation 
between the 880 /im source and the central K clump remains 
the same (3.2 ±0.2 kpc). 
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Fig. 2. — Lens modeling results. Major tickmarks are spaced at intervals of l^^ To ease comparisons, a cross is drawn at the center of each panel, a: Keck K 
image after subtracting Gl and G2. b: Best-fit K model convolved with the K-bmd PSF. Critical curves are in red and caustics are in blue. The box delineates 
the region covered by the source images (i.e., ej, and o). c: K model convolved with the SMA beam. It is clearly different from the SMA and EVLA images 
(/ & k), indicating differential magnification, d: K residual, e: Modeled intrinsic source morphology (i.e., without PSF; grey scale) vs. a direct inversion of the 
observed image (red contours). For comparison, the 880 /im (purple) and CO(l— ^'O) (green) sources are shown as color-filled ellipses. /: SMA 880 /im image. 
The grey ellipse shows the beam. Here and in i, contours are drawn at -2,-1,-1-1,-1-2, and +4a, where a is the r.m.s. noise (3 mJy beam~^). g: 880 /im model, h: 
Model convolved with the SMA beam, i: 880 /im residual, j: 880 /im source. The purple circle shows the FWHM of the source, k: EVLA CO(l— ^►O) image. 
Here and in n, contours are drawn at -1,-1-1,-1-2,-1-4, and -l-8cr, where a is the r.m.s. noise (27 /iJy beam~^). /: CO model, m: Model convolved with the EVLA 
beam, n: CO residual, o: CO source. The green ellipse shows the FWHMs of the source. 



3.4. CO(l^O) Source 

We use t he sa me technique to model the EVLA CO(l^O) 
map as in § 13.31 The lensed images are better resolved in the 
EVLA image than in the SMA image, so we use an elliptical 
Gaussian instead of a circular Gaussian for the source profile. 
The model has a total of six free parameters (x^y, FWHM, q, 
PA, and flux density). Again, we can constrain the astrometric 
offset between EVLA and Keck through lens modeling. Fig- 
ure [S? shows the minimum values at each offset position 
relative to the K image. To deal with the correlated noise, we 
scale the values from the residual map by the product of 
the FWHMs of the major and minor axes of the beam. The 
global best-fit, with reduced x^ of unity, is reached when we 
shift the EVLA image 0.''2 E of the K image. The bottom 
panels of Fig. [2] show the best-fit model. We estimate a CO 
magnification of /ico = 6.9 ± 1.6. Similar to the SMA image, 
we find a source-plane separation of 4.7 ± L6 kpc between the 
cold molecular gas and the stellar emission (i.e., the central K 
clump). The CO(l^O) is emitted from a more extended re- 
gion than the dust, but the two spatially overlap (Fig. [2^). The 
CO source has FWHM = 0'.'9 ± 0.''3 = 6.8 ± 2.3 kpc along the 
major axis, with an axis ratio of 0.8;^q 

The molecular gas disk is massive. The velocity-area- 
integrated CO brightness temperature of L(^q = [6.4 ± LO] x 
10^^ K km s~^ pc^ indicates a molecular gas reservoir 



of Mgas = [7.4 ±2.1] X 10^^ M© after lensing correction, 
assuming a conversion factor of aco = ^gas/^co = 0.8 
Mq (K km s~^ pc^)~\ which is commonly assumed for star- 
burst environments ( Solomon & Vanden Bout 2005). Note, 
however, that aco is uncertain by a factor of a few. 

The CO magnification factor determined from the lens 
model (6.9 ± L6) is in excellent agreement with that estimated 
from the CO luminosity-FWHM correlation. The observed 
Lqq and line width indicates a magnification factor of 7 ± 2, 
based on its deviation from the correlation established by un- 
lensed SMGs (Harris et al. 2012, Both well et al. in prep). 
This agreement demonstrates that strongly lensed SMGs may 
be effectively selected with CO spectroscopy in the future. 

4. SPECTRAL ENERGY DISTRIBUTIONS 

Useful physical parameters are encoded in the SEDs. The 
optical-to-NIR SED of HATLAS 12-00 is dominated by the 
foreground galaxies Gl and G2, from which we can derive the 
photometric redshift and the stellar population of the lensing 
galaxies. The far-IR and submillimeter regime is dominated 
by the lensed SMG, as evident in the SMA image, so the data 
can tell us the dust and star formation properties of the SMG. 

4.1. Lensing Galaxies 

Adopting the photometric redshift of 1.06, we model the 
nine-band photometry (u, g, r, /, z, Y, J, H, and K) 
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Fig. 3. — Constraining the astrometric offset between SMA, EVLA and Keck, a: The background image is the map of the best-fit models as a function 
of SMA-Keck offset. The image is displayed in logarithmic scale. Positive offsets indicate shifting the SMA image W or N relative to the Keck image. Iso-x^ 
contours are overlaid for 1, 2, and 3cr intervals. The white plus sign indicates the offset that yields the minimum value. The ellipse shows the Icr astrometry 
uncertainty determined from FIRST-SDSS cross-correlation. The cross indicates the systematic offset {-()'.' 01, -i-0f^08) between FIRST and SDSS within 1° 
radius of HATLAS 12-00. The overlapping area between the Icr contour of the x^ map and the ellipse gives the best estimate of the astrometric offset and its 
uncertainty. Zero offset is indicated by the white circle. b\ Same as a but for EVLA relative to Keck. 




Observed Wavelength (/xin) Observed Wavelength (/xnn) 

Fig. 4. — Modeling the SEDs. Top axes indicate wavelengths at the rest frame of the SMG (z = 3.26). a: Black data points are for the foreground lenses G1-I-G2, 
and the red data points are for the lensed SMG. The top black curve shows the best-fit BC03 stellar population synthesis model of G1-I-G2 at z = 1.06, using 
the nine data points below 3 fim. The bottom black curve shows the best-fit BC03 model for the SMG, along with the Icr range of acceptable models, b: The 
full SED of the SMG. The short dashed {green), dash-doted {purple), and long dashed {orange) curves are the best-fit SED templates of the "Cosmic Eyelash" 
flvison et al. 2010b), Arp 220, and Mrk 231, respectively. The Eyelash provides the best description of the overall SED among the three. The solid black and 
blue curves are the best-fit models with a single-temperature modified blackbody using the general and optically thin formula, respectively. The inset shows the 
1, 2, and 3cr contours in the T-^ plane for the general {black) and optically thin {blue) models. 



of G1+G2 with the stella r population synthesis n iodels of 
iBruzual & CharlotI (I2003L BC03). We assume a iChabrieri 
(l2QQ3h initial mass function (IMF), 'Calz etti et all (|1994) ex- 
tinction law, and exponentially declining star formation his- 
tory, with a range of e-folding times (r = 0.1 to 30 Gyr) and 
ages (0.01 to 12.5 Gyr). For each template, we fit for the 
stellar mass (Msteiiar) and extinction {E{B-V)). The best-fit 
model gives = 6.4 for dof = 7 (Fig.|4]2). The derived prop- 
erties of G1+G2 are listed in Table [3l The intrinsic extinction 



is small iE{B-V) = 0.04;^q;q4) and there is very little current 
star formation (SFR = 0.1;^q ^ Mq yr~^). The dust-absorbed 
UV/optical luminosity ([3^^^] >< ^qIO l^) is less than 0.15% 
of the total integrated IR luminosity before lensing correc- 
tion (Lg-iooo = 1.2 X 10^"^ L0). Therefore, Gl and G2 do not 
contribute significantly to the far-IR fluxes, in agreement with 
their absence in the SMA image. The stellar mass from SED 
modeling is ~80% of the total mass within the critical curves 
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TABLES 
Derived Physical Properties 



Object 


Quantity 


Value 


Unit 


G1+G2 


^stellar 

Age 
r 

E(B-V) 
SFRopt 


3.5+j-8 X IQii 
3±2 
0.3 ±0.2 
04+^-^^ 
1+0-4 


M© 
Gyr 
Gyr 
mag 
Mq yr-1 


SMG^ 


^stellar 

SFRopt 


[3.5 ±2.4] X 10^0 
1000 ±1000 


Mo 
Mo yr 1 


SMG^ 


^ust 

/3 

^dust 

Ao 

cr 

SFRiR 


62 ±3 
1.6±0.4 
[6.0±1.8] X 10^ 
210±40 
0.8 ±0.2 
[1.6±0.3] X 10^3 
1900 ±400 


K 

Mo 
/im 
kpc^ 
Lo 
Mo yr 1 


SMG'^ 


^dust 
^dust 


44±3 
1.3±0.2 
[8.4 ±2.6] X 10^ 


K 

Mo 


SMG^ 


L^Q [9.3 ±2.6] X IQio 
Lir/4o 172 ±58 
FWHMco 585 ±55 
Mgas [7.4 ±2.1] X 10^0 
Mdyn [3.2 ±1.3] X 10^1 

Mgas/Mbaryon 68 ± 17% 
Mgas /Mdyn 23 ±11% 


K km pc^ 
Lo/Kkm pc^ 
km 

Mo 
Mo 



Note. — Magnification-dependent parameters have been de- 
magnified and their errors include the magnification uncertainties. 
a - Stellar population synthesis modeling of the near-IR SED. b - 
General "optically thick" modified blackbody fit to the far-IR-to-mm 
SED. c - Optically thin modified blackbody fit. d - Parameters de- 
rived from CO(l— ^0) observations. 

from lens modeling (Table [3]) , implying that the galaxies are 
dominated by stellar mass within ~ 7 kpc. 

4.2. Lensed SMG 

In the / and ^-bands, we obtain the photometry of the 
SMG with an aperture contoured around the multiply-imaged 
features after subtracting the foreground lenses. We mea- 
sure 7^ = 21.2 ±0.1 and (J-K) = 2.1 ±0.2 in AB magni- 
tudes, consistent w i th the vedJ-K colors of unlensed SMGs 
(iFraver et al.l 12004 iDann erbauer et al."2004'). Careful mod- 
eling is required to extract SMG photometry from the WISE 
data, because the SMG is blended with the foreground galax- 
ies Gl through G4 (FWHM = 6-12'0. We model the WISE 
3.4 /im source with fo ur elliptical G aussians of the same 
shape. With GALFIT (Peng et all2010 ). we fix their positions 
to those determined from the Keck image, but we allow the 
Gaussian shape to vary. Then we measure the flux density 
of G1+G2 (~205 /iJy) decomposed from that of G3 and G4 
(-117 /iJy). Finally, the flux density of the SMG (-37 /iJy) 
is estimated from the excess of G1+G2 relat ive to their best- 
fit stellar population synthesis model (§ 14.11) . Unfortunately, 
we can not separate G3+G4 from G1+G2 in the longer wave- 
length channels of because of the inferior image quality, 
so we treat the excesses over the best-fit model of G1+G2 as 
upper limits for the SMG. The far-IR-to-millimeter SED is 
dominated by the lensed SMG, therefore no foreground sub- 
traction is necessary. Table [T] summarizes the photometry for 
HATLAS 12-00. 



We opt to model the rest-frame optical and far-IR emission 
separately, instead of fitting them together i n a self-consistent 
way with MAGPHYS (Ida Cunha et al.ll2008|) . because our lens 
model shows that they are emitted from physically distinctive 
regions, i.e., the dust that attenuates the optical emission has 
little to do with the starburst-heated dust that emits in the far- 
IR. 

It is difficult to constrain the stellar population with only 
three photometric detections (dof = 1) in rest- frame optical. 
However, we can limit the parameter space by excluding un- 
physical models, such as those that require negative extinc- 
tion corrections and those whose ages exceed the cosmic age 
at z = 3.2592 (1.9 Gyr). We use the same BC03 templates 
as in the previous section. The shaded region in Fig. ^ 
shows all of the permitted models with < dof + 1 . These 
models give a range of extinctions, stellar masses, and SFRs: 
E(B-V) < 0.94, Msteiiar = [3.5 ± 2.4] x 10^^ Mq, and SFR 
< 2000 Mq yr-i (TableEJ). The dust-absorbed UV/optical lu- 
minosity range from 0% to 140% of the observed Lg-iooo; but 
90% of the models have dust-absorbed luminosity less than 
50% of the observed Lg-iooo- 

We fit the far-IR SED with a single-temperature modified 
blackbody, 

^(^obs) = CT (l-e-^) 5(^rest, T) (1 +z) (3) 

where a is the total absorption cross section of dust par- 
ticles at the optically thick limit (i.e., the size of the 
dust-obscured region), B{y^T) the Planck function, r = 
(^rest/^o)^ = (Ao/Arest)^ the Optical depth, /i the tensing mag- 
nification factor, and di the luminosity distance. In the opti- 
cally thin limit (A ^ Ao), dust mass can be derived based on 
the knowledge of the opacity (absorption cross section per 
unit mass): 



/^^/(^rest) ^(i^rest, T) (1 +z) /i 

It is generally assumed that the opacity follows a power 
law, ndi^) oc , and has a normalization of i ^d = 0.07 ± 
0.02 m^ kg-^ at 850 (im (iDunne et al.1 120001: iJames et al.l 
120021) . Both the general "optically thick" {S^ oc (1 - 
e~^)By(T)) and the optically thin {Sy oc u^ByiT)) models pro- 
vide good fits to the observed SED (Fig.ll^?). For the general 
model, we use all of the nine detections between 100 /im and 
3 mm. The best-fit general model gives = 5.6 for dof = 5. 
For the optically thin model, we exclude the PACS 100 /im 
point, which is clearly on the Wien tail where small grains 
tend to dominate the emission. The best-fit optically thin 
model gives = 5.1 for dof = 5. The derived parameters 
are listed in Tabled 

The optically thick model yields dust properties simi- 
lar to those of the loca l Ultra-L uminous Infrared Galaxy 
(ULIRG) Arp 220 (RangwalaeLal. 2011), with the optical 
depth exceeds unity below —210 /im. The intrinsic 8 - 
1000 /im luminosity of Lg-iooo = [1.6 ±0.3] x 10^^ clas- 
sifies HATLAS 12-00 as a hyper- luminous infrared galaxy 
(Hy-LIRG). T he IR l uminosity impHes an SFR of 1900 ± 400 
Mq yr-i for a lChabri er (2003) IMF (Kennicutt 1998). Using 
the values of Lg-iooo and T = 62 ± 3 K in the Stefan-Boltzmann 
law we obtain a spherical source radius of 790 ±100 pc, 
which is three times larger than that of Arp 220 (230 pc) be- 
cause of the ten times greater luminosity. The source radius is 
comparable to that we derive from the optically thick model 



10 



Fu et al. 



(r = 510 ± 70 pc) and is consistent with the size we measure 
from modeUng the SMA image (FWHM = 1.2!i;^ kpc; §[33]). 
Therefore, the optically thick model is preferred. 

The radio luminosity from the observed 1.4 GHz flux 
density is Li.4ghz = 47rJ2^i.4GHz(l + z)""^ = [7.3 ± 3.4] x 
10^^ W Hz~^ for a radio spectral index of a = 0.7. As- 
suming the radio emission is magnified by the same fac- 
tor as the submillimeter emission, the SMG falls precisely 
on top of the radio-far-IR correlation of SMGs with ql = 
log(LiR/(4.52THz Li.4ghz)) = 2.1 ± 0.2 (see iKovacs et al.1 
I2OO6 ), suggesting no AGN contribution to the radio emission. 

We also do not see significant AGN contribution in the mid- 
IR. In Fig.gb, we fit the SEDs of the local ULIR Gs Mrk231 
Arp 22 0, and the z = 2.3 SMG "Cosmic Eyelash" (llvison et al.l 
|2010b) to the far-IR SED of H ATLAS 12-00. The upper 
limits He well below the AGN-dominated ULIRG Mrk 231 
but are more consistent with Arp 220 and the Eyelash. There- 
fore, we conclude that HATLAS 12-00 does not show any sig- 
nificant AGN contribution. 



5. DISCUSSION AND CONCLUSIONS 

We have presented high-resolution ^-band, 880 /im, and 
CO(l^O) observations and the near-IR-to-centimeter SED 
of a /f^r^c/z^/-selected strongly lensed SMG at z = 3.2592 
(HATLAS 12-00). The SMG shows distinctly different mor- 
phologies in the three images, suggesting differential magnifi- 
cation due to stratified morphologies. A joint strong lens mod- 
eling shows that the SMG is lensed by four galaxies at z ~ 1 
and the luminosity-weighted magnification factors are 16.7 ± 
0.8 in K, 7.6 ±1.5 at 880 /im, and 6.9 ± 1.6 at CO(l^O). In 
the source plane, the SMG consists of several stellar clumps 
extended over ^1.6 kpc with [3.5 ± 2.4] x 10^^ M© of stars, 
a compact (~1 kpc) starburst enshrouded by [6.0± 1.8] x 10^ 
of dust at ~60K, and an extended (~6 kpc) cold molecu- 
lar gas reservoir with [7.4 ±2.1] x 10^^ M© of gas. The 
starburst and its gas reservoir are located ~4 kpc from the 
stars. Similar separations between optical and submillimeter 
emission have been observed in unlensed z> 2 SM Gs (e.g., 
iTacconi et al.ll2008l: ICmUTetan 2010: R iechers et al. 2010). 
However, these previous results could also be attributed to as- 
trometry offsets across different facilities. Here the physical 
separations amongst stars, dust, and gas are less ambiguous 
because of the clear wavelength-dependent morphologies in 
the image plane: fortuitously, the stars and dust/gas strad- 
dle across the caustic, so the less obscured /^-band region is 
quadruply imaged while the heavily obscured starburst and its 
gas reservoir are doubly imaged. Because of the ~4 kpc sep- 
aration between the stars and the gas-rich starburst and their 
similar masses, it is tempting to suggest that the SMG is in 
the process of a major merger, which presumably is driving 
the starburst activ ity in z > 2 SMGs (e.g., Tacconi et al. 2008; 
lEngel et ani2010l) . However, spatial separation of this scale 
can also be explained by differential dust obscuration in a sin- 
gle galaxy, as has been proposed for other high-redshift SMGs 
where spatial offsets have been observed between rest-frame 
UV and submillimeter (e.g., GN20 and AzTEC 3; CarilH et al] 
[MlOHRiechers eta l. 2010). 

How does the massive gaseous disk compare with the 
disks in other SMGs? We estimate an dynamical mass of 
[3.2±1.3] X 1O^^M0 and a gas fraction of /gas =Mgas/Mdyn = 
23 ± 1 1% for the C0(1— ^0) disk using the "isotropic virial es- 



timator" (e.g.. ITacconi et al.ll2008"h : 

Mdyn = 2.8 X lO^AV^HM^HWHM M©, (5) 

where AVfwhm is the CO line FWHM in km s~^ and rnwHM 
is half of the FWHM size of the disk in kpc. Combined 
with the FWHM disk radius, we further estimate a gas sur- 
face density of S = 510 ± 370 M© pc~^. Both /gas and S are 
similar to those of the extended CO(l^ O) disks in the two 
z-^3A SMGs in Riechers e t al.l (1201 Ibl) . But both values 
are significantly smaller than those of the kinematically re- 
solved CO(6^5) disk of t he "Cosmic Eyelash" (/ gas ^ 70%, 
S - 3000 ± 500 Mo pc-2; iSwinbank etani2011h . The dis- 
crepancies illustrate the limitations of these wildly used but 
crude estimators and/or that high excitation CO lines probe 
more compact and denser regions in a disk. Higher resolution 
observations are clearly needed to resolve this issue. Without 
spatially resolved gas kinemat ics, we refrain from estimating 
the disk stability parameter of Toomrd (119641) . 

HATLAS 12-00 is a gas-rich, initial starburst system sim- 
ilar to unlensed SMGs and local ULIRGs. Its intrinsic IR 
luminosity well exceeds 10^^ L© (i.e., Hy-LIRG), implying 
an enormous rate of star formation (1900 ±400 M© yr~^). 
Although the molecular gas reservoir is massive and it con- 
stitutes 68 ± 17% of the visible baryonic mass (Mgas +Mdust + 
^stellar) and 23 ± 11% of the dynamical mass, it will exhaust 
in just 39 ± 14 (aco/0.8) Myr at the current SFR (assum- 
ing no gas accretion). The star formation timescale, tsf = 
^steiiar/SFR = 18 ± 13 Myr, is only ~1% of the cosmic age at 
z = 3.2592 (Tcosmic = 1.9 Gyr), suggesting that HATLAS12-00 
is an initial starburst system with maturity /i = tsf/ Tcosmic <C 1 
(iScoville et al.ll2007h . The specific SFR, sSFR = SFR/Msteiiar 
= 54 ± 38 Gyr~\ is consistent with the average z> 2 SMGs, 
but it is an order of magnitude higher than the median value 
of the star-forming main sequence of Lyman break galaxies at 
the same epoch (Daddi et al. 2009). The dust mass of [6.0 ± 
1.8] X 10^ M(7) is similar to the average dust mass of unlensed 
SMGs (e.g., Michalowski et al. 2010), and the gas-to-dust ra- 
tio, Mgas/Mdust = 120 ±50 (/^850Mm/0.07 m2 kg-') (aco/0.8), 
is comparable to that of the Milky Way. Assuming that the 
dust emission is indicative of the size of the starburst and a 
starburst disk radius of ro = 1 kpc, the star formation surface 
density of 600 ±120 Mq yr~^ kpc~^ approaches the Ed- 
dington limit of radiation pressure supported starburst disks, 
similar to local ULIRGs (t^ - 10^ M© yr'^ kpc-^i lScovilE 
2003; Thompson et al. 2005). 

In conclusion, HATLAS 12-00 is a bona fide SMG with an 
intrinsic submillimeter flux density of 5*880 = 9.2 ± 2.2 mJy. 
The starburst disk, where most of the molecular gas and dust 
reside, is spatially separated from the less obscured stellar 
population by ~4 kpc, suggesting either a major merger or 
differentiated dust obscuration. The ~1 kpc radius starburst 
disk is presumably supported in large by radiation pressure on 
the dust grains. Its physical properties, such as molecular gas 
mass, stellar mass, gas-to-dust ratio, gas fraction, SFR, star 
formation efficiency, and radio-to-far-IR luminosity ratio, are 
all very simi lar to unlense d z>2 S MGs ("Hainline et al."201 1"; 
Wardlow et a l. 2011; Michalowski et al. 2010; Kovacs et al.J 
20Q6h . The lensing boost of the effective angular resolution 
and sensitivity has allowed us to examine in unprecedented 
details the properties of a typical starburst galaxy when the 
Universe is only 1/7 of its current age. HATLAS 12-00 pro- 
vides a prelude to a golden age of SMG research, as Herschel 
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is unveiling hundreds o f strongly lensed SMGs bef ore the 
mission completes (e.g.. iGonzalez-Nuevo et al.ll2Q12L Ward- 
low et al. in prep). 
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